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Dynamical tunneling has been observed in atom optics experiments by 

in 

£S| . two groups. We show that the experimental results are extremely well de- 

o , 

scribed by time-periodic Hamiltonians with momentum quantized in units of 
. the atomic recoil. The observed tunneling has a well defined period when 



X 



only two Floquet states dominate the dynamics. Beat frequencies are ob- 



served when three Floquet states dominate. We find frequencies which match 

a ; 

those observed in both experiments. The dynamical origin of the dominant 



Floquet states is identified. 
PACS numbers : 03.65.Xp, 03.75.Be, 05.45.Mt, 42.50.Vk 

Atom optics experiments recently have been used to investigate the effect of underlying 
classical chaos on quantum dynamics. The experiments we focus on in this paper have 
demonstated the existence of dynamic tunneling in momentum space in regimes where the 
underlying classical phase space contains a mixture of chaotic and regular orbits. We will 
show that we can accurately reproduce the dominant tunneling frequencies observed in these 
two very different experiments using Floquet analysis of the quantum dynamics. 

Typically, cold sodium or cesium atoms are allowed to interact with laser beams which 
are detuned away from resonance with two atomic energy levels which have energy spacing, 



hujo- Two counterpropagating laser beams create a periodically modulated standing wave 
of light which stimulates absorption and then emission of a photon. This results in a net 
atomic recoil of 2hki, where kt, = ujl/c is the wave vector of the laser beams and H is 
Planck's constant. When the laser detuning, 5l = c^o — is large, this process dominates 
the dynamics. 

A theoretical model which describes the atomic dynamics in such systems was developed 
by Graham, Schlautmann, and Zoller [1]. Recently, two groups, Steck, Oskay and Raizen 
[2,3] in Texas and Hensinger et.al. at NIST [4], have performed independent experiments 
in which dynamic tunneling has been observed. In this letter, we explore the accuracy of 
the models used to analyse these experiments, and the dynamical origin of the tunneling 
observed in each experiment. We first discuss the Texas experiment and then the NIST 
experiment. 

In the Texas experiment [2,3], the dynamics of non-interacting cold cesium atoms, in 
an amplitude modulated standing wave of light, was measured. The atomic center-of-mass 
Hamiltonian (in S.I. units) used to model dynamics of the cesium atoms is 

^ = £ " 2Vo cos2 (ir) cos {2kL$) ' (1) 

where p, x, and m are the momentum, position, and mass, respectively, of a cesium atom, 
u m = is the modulation frequency, and V Q = H ^J^ X is the ac Stark shift amplitude, where 
^max = —2E d/h is the Rabi frequency, E is the electric field strength, and d is the dipole 
moment of cesium [5]. 

In the experiments, the initial state is well localized at discrete momentum states sep- 
arated by 2Tiki. This quantization of the momentum occurs naturally in the experiment 
due to the presence of counter propagating laser beams which cause two-photon transitions. 
Therefore, in our theoretical analysis, we perform a scaling which explicitly quantizes the 
momentum in units of 2hki- As we will see later, this allows us to use Floquet theory rather 
than Floquet-Bloch theory deals with a continuum of monentum states [6]. Let if> = 2kix, 
p = 2nhkL, uj r = %k 2 L /2m, u = uj m /Auj r , t' = 4u r t, and H th = mH /2k 2 L h 2 , to obtain, 
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cos(0) H — cos(0 — out') H — cos(0 + ujt') 
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where a = 8u r T 2 V /h. All quantities are dimensionless and h is the dimensionless momen- 
tum operator with eigenstates, \n), and integer eigenvalues, — oo<n<oo. (Note that the 
experimental papers [2,3] perform the following scaling, = 2kix, r = u m t/2n = t/T, 
p = AirkLp/mUm, H exp = 16n 2 k\H / mw^, and Hamiltonian takes the form, H exp = 
y — 2acos 2 (7rr)cos(0).) This system has three primary resonances centered at (n — 0, <p = 0) 
and (n = ±cj/2, = 0). For small values of a (a < 1.5), the primary resonances have 
pendulum-like structure, and the resonance at n = has half-width An = \J f^r, while the 
resonances at n± = ±cj/2 have half-width An± = An /^/2 [7]. The primary resonance at 
n = bifurcates at a~7.0. The two outer primaries remain visible until a~13.0 when they 
disappear. 

The classical motion is obtained from Hamilton's equations, n = — |^ and <j) = ^§^- 
In the Texas experiment, u r = 1.30xl0 4 rad/ s and T = 2n/uj m = 20/is, which, for small 
a, gives a location of n± = ±3.0 for the outer primary resonances. For the field strength 
a = 9.7, used in the Texas experiment, the pendulum approximation for the half-widths of 
the two outer primary resonances gives An± = 2.1, while the half-width of the central island 
is An = 3.0. Thus the Texas experiment, for a = 9.7, is in the strong field regime, where 
the primary resonances have overlapped and considerable chaos is expected [7] . A surface of 
section of the classical phase space for a = 9.7 is shown in Figure (l.a). The central primary 
resonance has bifurcated and is largely destroyed, and the outer primary resonances have 
been reduced significantly in size and are centered at momentum values n = ±4.2. Note 
also that the chaotic region lies in the interval — 5<n< ±5, indicating that eleven quantized 
momentum states determine the dynamics in the chaotic region. 

The Texas experiment used atoms prepared initially with a narrow momentum distribu- 
tion peaked at n = 4.2 (on the upper island). To numerically simulate this initial condition, 
we solved the Schrodinger equation, i ^Jf, - = H t h\^(t')), using momentum states, |n), as a 
basis. A coherent state, 



(n|*(0)>= (n \<p o n ) = — exp 
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centered at (n = n D , = O ) is used as the initial state, with cr = 1.2 which was used in 
the experiment. In the momentum basis, the Schrodinger equation reduces to a system 
of coupled first order differential equations for the amplitudes, {n\^f(t)). This system was 
truncated, and 81 equations for the states (n\^(t)) with — 40<n<40 were kept. The time 
variation of the average momentum, (n), is shown in Figures (2. a), (2.b), and (2.c) for 
a = 8.0, 9.7, and 13.0, respectively. In all cases, the initial state is (n = 4.2, <p = 0). For all 
three plots, the average momentum oscillates between the outer primary resonances. The 
plot for a = 8.0 (a = 9.7) has two dominant frequencies, f\ = 1.95kHz and / 2 = 2.73kHz 
(/i = 2.39kHz and f 2 = 2.88kHz), giving rise to a beating effect. The beating effect at 
a = 9.7 was observed in the Texas experiment [3], but the experimental error bars were 
too great to resolve it at a — 8.0. The plot for a = 13.0 shows one dominant frequency, 
/ = 1.56kHz. 

It is useful to examine these results using Floquet theory [7]. Because the Hamiltonian, 
H thl has time periodic coefficients, the Schrodinger equation has Floquet solutions of the 
form (n\^(t)} = e~ tnjt {n\xj(t)) where f^- is the j th Floquet eigenphase and \Xj(t)) is the j th 
Floquet eigenstate and is periodic in time, \xj(t)) — \Xj(t+T)) [7]. The Floquet eigenphases, 
Qj, are conserved quantities, and the eigenstates form a complete orthonormal basis which 
can be used to analyze the dynamics. The states, |Xj(0)), are eigenfunctions of the Floquet 
matrix, U(T), and the phase functions, e ~ lf ^ T , are its eigenvalues. The Floquet matrix is 
computed by by taking a momentum eigenstate as the initial state and evolving it for one 
period, T, using the Schrodinger equation. The resulting vector (in the momentum basis) 
is a column of the Floquet matrix. 

The overlap probabilities, P 7 = |(xj(O)|0 o n o (O))| 2 , give the contribution of each Floquet 
state to the dynamics. The probability to find the system in momentum state, \n), at time 
t, can be written [7] 
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I (n \<p n {t))\ 2 = ]T]T ex P H ( n i " n <) *) ( n IXi (*)) (X< (0 \n) 

i 3 

x(x,(O)|0 o n o (O))(0 o n o (O)| X ,(O)), (4) 

with time t in seconds and Qj/27r in Hertz. The oscillation frequencies, f exp , observed in 
the experiments can be equated to differences between Floquet eigenphases. The frequency 
differences, f exp = (Sl^- — O^) /27r, for Floquet eigenstates with overlap probability, PjPj>0.04, 
are plotted in Figure (3) for the range of parameters shown in the Texas experiment [3]. Each 
curve is a plot of the frequency difference between two Floquet states as a function of a. At 
values of a where there are multiple curves there are more than two dominant frequencies. 
The Texas experiment was able to resolve the dominant frequencies, f exp < 3kHz, in the 
interval between a~8.7 and a~10.3. Our analysis exactly reproduces those experimental 
results. In the amplitude range, a^7.6 to awll.6, we find that two frequencies dominate 
and give rise to the beats seen in Figures (2. a) and (2.b). In the Texas data [3], large error 
bars occur in the regions a = <7.0 and a>13.7. This may be due to the rapid change in 
the dominant frequencies in those regions. A fundamental change in the dynamics occurs 
for a > 14, where a different set of Floquet states begins to dominate the dynamics. 

Only eleven Floquet states have support on momentum in the region, n = — 5 to n = 5, 
and determine the dynamics in the chaotic region. In Figures (4.a)-(4.d), we show Husimi 
plots for the Floquet states which, for a = 9.7, have the largest overlap probability, three 
of which dominate the dynamics. The dark regions of these plots show the region of the 
classical phase space where the probability of finding the cesium atoms is largest. The 
eigenphase differences, (Q^ — ^4a)/2vr = 2.89/cHz and (fi^ — Q4 C )/2n = 2.40/cHz correspond 
to the two dominant oscillation frequencies observed by the Texas experiment at a — 9.7. 
The state in Figure (4.d) has the fourth highest overlap probability, P d = 0.045, but it lies 
in the chaotic sea. The state in Fig. (4d) and others not shown contribute to the fine scale 
structure in these curves. 

Let us now consider the NIST experiment [4], which used a Bose-Einstein condensate of 
sodium atoms to observe dynamic tunneling. Formation of a condensate with the sodium 
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atoms yields a narrower distribution of initial momenta than the Texas experiment. The 
Hamiltonian used to describe the experiment can be written in the form 

,~,2, 



H th = r? + — - 



1 + 2z/e cos(a>t') — cos <f> — z/ecos(0 — out') — z/ecos(0 + Cot') 



(5) 



where v = ±1. When v — +1 {y — —1), Eq. (5), with starting time t' = 0, 
reproduces the dynamics of the NIST experiment which has starting time, f = T/4 
(f = 3T/4). This Hamiltonian again discretizes the momentum in units of 2hki, which 
reflects the quantization of momentum due to the two photon transitions. (The Hamilto- 



nian, H exp = ^- + 2k (1 + 2esin(u; m r)) sin (|), used in the experimental paper is obtained 
by setting 3 = and H exv = 8fc 2 L - ? l H th .) 

For small amplitudes, k and ne, the NIST Hamiltonians have three primary resonances. 
For v — — 1 they are located at (n — 0, = 0) and (n = = ±7r), while for v — +1 they 
are located at (n = 0, = 0) and (n = ±|, = 0). They have half- widths, An = \Zuj 2 k 
and An± = \/uj 2 k€ [7]. 

A strobe plot of the classical phase space for the Hamiltonian in Eq. (5) with v — — 1 
and experimental parameters oj m /2n = 250kHz, Co = 2.5, k = 1.66 and e = 0.29 is shown in 
Figure (l.b). Seven Floquet states determine the dynamics in the chaotic region between 
n = —3 and n — 3. For the parameters used in the experiment, the pendulum approximation 
predicts the primary resonances to lie at n = and n = ±1.25, and have half-widths, 
An = 3.2 and An± = 1.7. We find that the primary resonances are totally destroyed at 
K&0.2, and then new resonances, which resemble the primaries, reappear and disappear 
repeatedly as k is increased. For k = 1.66 and e = 0.29, a large resonance exists at 
(n = 0, = 0) and three small pairs of higher order resonances exist at (n^±1.5, = ±7r), 
(nw±3.0, = 0) and (ra«±2.0, = 0). 

In Figure (5) we show the time evolution of the momentum expectation value for two 
different initial conditions for the v — — 1 Hamiltonian at parameter values, k = 1.66, 
e = 0.29, Co = 2.5 and Co m /2n = 250kHz. Figure (5. a), with (n = 1.6, O = 0), shows a 
somewhat noisy oscillation with a dominant frequency 24.9kHz (10.0 modulation periods), 
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which is in good agreement with the experimental result. Figure (5.b) shows the case with 
(n = 3.0, O = 0). A clean oscillation with frequency 18.3kHz (13.7 modulation periods) 
occurs. This oscillation was not observed in the experiment, but we expect it would show 
up in a power spectrum of the experimental data. 

We now consider a Floquet analysis for both Hamiltonians, v — ±1. The Floquet eigen- 
phases for v = ±1 are identical, but the Floquet eigenstates associated with each eigenphase 
are different for the two Hamiltonians. Let us first consider the v — — 1 Hamiltonian with 
parameters, k = 1.66, e = 0.29, Co = 2.5 and Lu m /2n = 250kHz. In Figures (6. a) and 
(6.b), we show the two Floquet states which dominate the dynamics for initial condition, 
(n Q = 1.6, 4> Q = 0). They have a frequency difference, (f2 66 — f2 6a )/27r = 25.0kHz. Their 
frequency difference accounts for the oscillation of 10 modulation periods reported in [4]. 
These Floquet states are not even-odd pairs as suggested in [4], and they both lie in the 
chaotic sea. If the effective Planck's constant for this experiment were smaller, more Floquet 
states would be supported by the chaotic region and we would not expect to find this simple 
oscillation [8] for this initial condition. 

If we take initial condition, {n = 3.0, <p = 0) for v = —1, we obtain the oscillation 
shown in Figure (5.b). This oscillation results from the even-odd Floquet pair shown in 
Figures (6.c) and (6.d). Figure (6.c) (Fig. (6.d)) is even (odd) under the transformation 
p^ — p. They have a frequency difference, (f2 6c — f2 M )/27r = 18.3kHz. This oscillation 
appears to result from states sitting the outer-most nonlinear resonance. 

We finally consider the v — +1 Hamiltonian with parameters, k = 1.66, e = 0.29, 
oj = 2.5 and Cj m /2n = 250kHz. We find that the 25.0kHz (10 modulation periods) oscillation 
dominates those initial momentum states which are centered at = and lie in the interval 
n = 1.7 to n = 3.0. These oscillations appear to result from the two Floquet states which lie 
in the chaotic sea. If we change the parameters to k = 1.82 and e = 0.30 and the modulation 
frequency to uj m = 222kHz, the dominant frequency for initial state, {n = 2.0, <p = 0), is 
36.8kHz (6.03 modulation periods), which is in agreement with the NIST experiment. 

In conclusion, the model Hamiltonians, with momentum quantized in units of 2hki, 

7 



give extremely good predictions of the experimental results. Because of the momentum 
quantization imposed by the dynamics of the experiment, we found that it was advantageous 
to use Floquet theory rather than Floquet-Bloch theory to analyse the experiment. In fact, 
our results are so good that these models might be used to help calibrate future experiments. 
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FIGURES 

FIG. 1. Classical strobe plots: (a) The Texas experiment with u = 6.0 and a = 9.7. (b) NIST 
experiment with to = 2.5, a = 1.66 and e = 0.29. 

FIG. 2. Evolution of average momentum, (n) (in dimensionless units) for the Texas experiment 
for uj = 6.0: (a) a = 8.0; (b) a = 9.7; and (c) a = 13.0. 

FIG. 3. Oscillation frequencies, AQ = (Qj — fij), calculated from the Floquet eigenphase dif- 
ferences for varying dimensionless field strengths, a. A threshold of PjPj>0.04 overlap probability 
was used to select the dominant frequencies. The three values shown at a = 9.7 correspond to 
(fLkj - 4b )/27r, (f2 4a - f2 4c )/27r, and (J2 46 - 4c )/27r 

FIG. 4. Husimi plots of Floquet eigenstates for the Texas experiment for u = 6.0 and a = 9.7. 
(a) Floquet eigenphase 4a /27r = 16.9kHz and an overlap probability, P 4a = 0.416. State (b) 
Floquet eigenphase r2 4 ft/27r = 19.7kHz and an overlap probability, P 4 & = 0.224. (c) Floquet 
eigenphase 4c /27r = 14.5kHz and an overlap probability, P 4c = 0.20. (d) Floquet eigenphase 
^4<i/2vr = 18, 389Hz and an overlap probability, = 0.045. 

FIG. 5. Evolution of momentum expectation value, (n) (in dimensionless units) for the NIST 
experiment for k = 1.66, e = 0.29, u = 2.5 and u) m /27r = 250Hz: (a) n Q = 1.6, and <p Q = 0; (b) 
n = 3.0, and <f) = 0. 

FIG. 6. Husimi plots of Floquet eigenstates for the NIST experiment with k = 1.66, e = 0.29, 
Co = 2.5 and uj m /2ir = 250Hz. (a) Floquet eigenphase Q,Q a /2TT = 49.0kHz and overlap probability 
PQa = 0.380. State (b) Floquet eigenphase Q,^/2tt = 73.9kHz and overlap probability Pqj, = 0.306. 
(c) Floquet eigenphase Qq c /2tt = 15.3kHz and overlap probability Pq c = 0.427. (d) Floquet 
eigenphase Qed/2TT = 33.5kHz and overlap probability P§d = 0.421. 
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